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It is shown that discontinuity, inherit to a square-well potential which is often used for mimicking 
the Ceo cage potential, results neither in qualitative nor, what is even more important, quantitative 
artifacts in problems of endohedral AQCeo atoms, such as their photoionization. Moreover, it is 
demonstrated that, upon mimicking a square- well potential by a potential with diffuse boarders, cal- 
culated photoionization spectra are largely insensitive to the degree rj of diffuseness of the potential 
boarders, in a reasonably broad range of t^'s. The proof is based on results of comparison between 
calculated data for HOCeo and XeQCeo photoionization cross sections and photoelectron angular 
asymmetries obtained by modeling the Ceo cage by confining square-well and diffuse potentials. 
Thus, it matters little which of these potentials to use in A@Ceo problems. The diffuse potential is 
modeled by a combination of two Woods-Saxon potentials. Calculated data for the photoionization 
of Xe@Ceo are separately detailed as well, since they are important stand-alone qualities. 

PACS numbers: PACS: 32.80.Fb, 32.30.-r, 31.15.V- 



I. INTRODUCTION 

Photoionization of atoms A encapsulated inside the 
hollow interior of the Ceo fuUerene cage, labeled as 
A@Cqo and referred to as endohedral or confined atoms, 
has been an ad hoc topic of numerous theoretical (see 
review papers P, (and references therein in addition 
to other references in the present paper) and experimen- 
tal [33 (and references therein) studies in recent years. 
An ab initio, based on first principles theory of A@Cgq 
photoionization has not been developed yet, in view of 
a formidable complexity of the problem. Perhaps, the 
most sophisticated theory here is represented by a time- 
dependent local density approximation (TDLDA), see, 
e.g., Ref. Q and references therein. However, TDLDA 
theory has obvious drawbacks in view of an unsatisfac- 
tory agreement with the existing experimental data on 
the Xe@C^Q dd photoionization [31]. Meanwhile, many 
important insights into the problem can be, and have 
been, unraveled on the basis of simpler empirical models 
based on the modeling of Cgo confinement by pseudo- 
potentials, such as a (5-function-like potential [5|, |6| (and 
references therein) or square- well-potential, /7swp(?')- 
The latter is defined as 



C^swp(»') 



-C/o, a Ro<r <Ro 
otherwise. 



A 



(1) 



Here, Rq is the inner radius of Ceo, ^ is the thickness of 
the Ceo wall, and Uq is the potential depth. The square- 
well-potential modeling of Cgo has become quite popular 
among various researchers and been used on numerous 
occasions in the field of Cen and A@Cqo theoretical stud- 
ies, see, e.g., Refs. [H, 0, □tIq'I and references therein. 

The potential Uswpif), however, is discontinuous at 
its boarders. A possible emergence of qualitative and, 
especially, quantitative artifacts, associated with said 
discontinuity, in calculated photoionization cross sec- 
tions of, and photoelectron angular- asymmetries from, 
A@Ceo atoms has not been detailed in literature ex- 



cept for a couple of brief remarks dropped on the sub- 
ject in passing in Refs. [l^, [2l|. In particular, it is be- 
cause of the lack of such knowledge that some authors, 
Ref. 



e.g. 



22| 



have recently proposed to discard 
the square-well-potential modeling of the Ceo cage on 
grounds of discontinuity of the potential, and replace 
it by a Gaussian-function-like model potential. This 
proposition, we argue, is wrong. The proposed poten- 
tial continuously changes strongly everywhere inside the 
Ceo cage, thereby having no compact boarders. This 
contradicts a recent Fourier imaging study of the ex- 
perimental Ceo photoionization cross section [23j . From 
the study, the fact that the Ceo has well-defined, sharp 
boarders is obvious. In short, following the logic line of 
Ref. [13, , photoionization of an atomic cluster occurs 
with the greatest probability where the cluster's poten- 
tial changes sharply. This is obvious from the acceler- 
ation form gage for a dipole photoionization amplitude, 
D, namely, D oc {ip f\ArV (r)\tpi) ■ On grounds of this, 
the Fourier imaging study [231 of the experimental Ceo 
photoionization cross section [3l| has proven that spheri- 
cal boundaries of Ceo are rather sharp. Hence, the Ceo 
potential should be reminiscent of a square-well poten- 
tial J7swp(^)j to a good approximation, rather than a 
soft Gaussian-function-like potential peaked at the Ceo 
mid-wall. 

Surely, a Ceo confining potential with diffuse (but well 
defined) boarders is much more realistic than a square- 
well-potential with infinitely sharp edges. Correspond- 
ingly, replacement of the latter with the former in A@Ceo 
problems would be, of course, worthwhile, for the sake 
of consistency. With this in mind, it is not at all clear 
apriori to what degree this may qualitatively and quan- 
titatively affect/alter results of numerous previous cal- 
culational predictions on A@Geo photoionization made 
on the basis of the square-well-potential concept, and 
whether such replacement will improve agreement with 
experiment or make it worse. The need for clarifying this 
issue for a deeper insight into the structure and spec- 



tra of AQCgo atoms is additionally signified by the fact 
that the square- well-potential model of A@Cqo atoms has 
been used abundantly over the years which has resulted 
in a large array of predicted data and phenomena that 
might need to be re-studied/re-calculated with an eye on 
more realistic diffuse potential boarders. 

It is the aim of this paper to dot i's and cross Vs con- 
cerning possible consequences of replacement of the con- 
fining square- well-potential Us-wpir) by a diffuse poten- 
tial Uj^pir) in A@Cqo problems. To achieve this end, we 
calculate the photoionization cross sections of, and pho- 
toelectron angular distribution from, inner and valence 
subshells of AQCgo atoms by utilizing both Uswp{t) and 
Uupir) [which, nevertheless, is reminiscent of Uswp{r), 
see next section] as model pseudo-potentials that mimic 
the Ceo potential. It is proven in the present work that 
discontinuity of the C/swp('') potential in question does 
not matter at all, in photoionization problems. It leads 
neither to qualitative nor quantitative artifacts in cal- 
culated photoionization characteristics of AQCgo atoms 
compared to calculated data obtained with the use of 
the diffuse potential UDp{r). Moreover, it is found that 
the degree of diffuseness of the C/dp('') potential matters 
surprisingly little - a fraction-to-nothing - in a relatively 
broad range of its values. 

Atomic units (a.u.) are used throughout this paper. 



II. THEORY, RESULTS, AND DISCUSSION 
A. Review of theory 

In this work, to model the Cgo cage by a square-well- 
potential we use Uswp{t) given by Eq. ([T]). 

To mimic the Ceo cage by a diffuse potential Upip{r), 
which is somewhat reminiscent of Uswpif), we form the 
latter by the appropriately defined combination of two 
Woods-Saxon potentials: 



C/dpW 



2Uo 



l + exp(^) 

2Un 



r<i?o + iA 



(2) 



Here, ry is the diffuseness parameter, and i?o, Uq, and A 
are the square-well-potential parameters. 

Concerning photoionization of AQCeo atoms, we fo- 
cus on the photoionization cross section crn£(w) of a 
n£-subshell of A@Ceo as well as dipole photoelectron 
angular-asymmetry parameter Pneiui) which governs the 
pattern of a n£-photoelectron angular distribution, see, 
e.g., Refs. [13, [23 fo'" details. Expressions for ani{u>) and 
fUnii^) for A@Cqo are the same as for free atoms pre- 
sented in the above references. We write, 



3(2^4-1) 



Lu[£\D,_i\^ + {£ + l)\Di+i\' 



(3) 



and we recast l3ni{uj) as 

e{£ - i)p^ - 6e{e + i)pcos^ + {i + 1){£ + 2) 



where 



P 



(2£+l)(p2^ + ^+l) 
\Di-i\ 



-1- 



,(4) 



(5) 



In the above equations, uj is the photon energy, Nne is 
the number of electrons in a n£ subshell, a is the fine- 
structure constant, Di±i are radial dipole photoioniza- 
tion amplitudes, and Si±i are phase shifts of Di±i. Note 
that the quantities ct„^, /3„£, D^±i, p, 6i±i, and <& all de- 
pend upon photon energy w; the explicit dependence is 
omitted in the above equations for reasons of simplicity. 
In a one-electron approximation, Di±i is defined as 



{r)rPni {r)dr, 



(6) 



where Pni (r) and P^^ (r) are one-electron radial wavefunc- 
tions of the bound and continuous states, respectively. 
In the present work, these wavefunctions and energies 
of a discrete spectrum are the solutions of a ra- 
dial Hartree-Fock (HF) equation accounting for a J/ceo 
confinement: 



[H^ + Uc,Ar)]Pnt{r) 



(7) 



Here, is a radial part of the HF Hamiltonian which 
is identical to that for a free atom, see, e.g., Ref. [23| for 
the latter, and Uceo (^) is either the square- well-potential 
Uswp{r), Eq. or diffuse potential Uup{r), Eq. In 
a particular case of an endohedral single-electron hydro- 
gen atom, HQCeo, Eq. ([7]) reduces to a radial Schrodinger 
equation in the presence of the Coo confinement: 



1 d^Pnje) 

2 dr^ 



-1 I ije + i) 

r 2r'^ 



= En(e)lPn{c)l{r)- (8) 



To calculate Di±i beyond a one-electron HF approxi- 
mation, i.e., to account for impacts of electron correlation 
on the photoionization process, we utilize a random phase 
approximation with exchange (RPAE) [2^ ; the reader is 
referred to Ref. [2^] for details of the RPAE theory. 



B. Results and discussion 

In this work, when approximating the Coo potential by 
a square-well or diffuse potential, we use the proposed in 
Ref. dil values of Uq = -0.422, A = 1.25, and Ra = 6.01 
a.u. rather than the earlier used Uq = —0.302, A = 1.9, 
and i?o — 5.89 a.u., see, e.g., Refs. 111. [l3 . This is 
because the new parameters were shown \2(y\ to result 
in a much better match to experimental photoionization 
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FIG. 1: (Color online) The square-well-potential f/swp(''), 
Eq. (O, with (7o = -0.422, A = 1.25, and i?o = 6.01 a.u. 
(dotted line), as well as the diffuse potential I7dp(''), Eq. ((2|, 
with the diffuseness parameter r\ = 0.01 (solid line) and r\ — 
0.1 (dashed line). See text for details. 

spectra of endohedral atoms (of Xe@C^Q Q , in particu- 
lar). 

Calculated data for the thus defined potentials 
C^swp('') and Uupir) are compared with each other in 
Fig.Ufor two different values of the diffuseness parameter 
77 of Uup{r), namely, 77 = 0.01 and 0.1 . One can see from 
Fig. [1] that Uup{r) with 77 = 0.01 matches the squarc- 
well-potcntial practically precisely. On the contrary, the 
substitution of 77 = 0.1 in Uup{r) results both in rela- 
tively strongly diffuse boarders and bottom of Uupir), 
so that the latter noticeably deviates from the square- 
well-potential. 

In the following, we detail results of comparison be- 
tween ane{uj) and Pnii'-o) calculated in the presence of 
the above determined J7swp(^) and Up,p{r) confinement. 

As the first step, in order to avoid various possible com- 
plications associated with electron correlation in multi- 
electron atoms, let us discuss a "clean" case - the pho- 
toionization of the Is ground-state of the endohedral 
hydrogen atom, HOCeo. Corresponding results are de- 
picted in Fig. [5] First, we note that the photoionization 
cross section o'is(w) of HQCeo calculated with the use 
of the square- well-potential C/swp (f ) and diffuse poten- 
tial UY)p{r) with a small diffuseness parameter 77 — 0.01 
are undistinguishable from each. Hence, discontinuity of 
the square-well-potential does not matter at all in this 
problem. Second, what comes as a big surprise, how- 
ever, is that calculated data obtained with the use of a 
much greater diffuseness parameter, 77 = 0.1, for which 
C^Dp(f) deviates from a square- well-potential consider- 
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FIG. 2: (Color online) Hydrogen Is photoionization cross sec- 
tion of HQCgo calculated with the use of the square-well- 
potential C/swp(7') (open circles), diffuse potential Uwir-) 
with ri — 0.01 (solid line) and r; = 0.1 (dashed line), and 
of free hydrogen, as marked. 



ably (see Fig. [TJ, are practically the same as calculated 
results obtained with 77 = 0.01 or with the use of the 
square-well-potential. This is the initial indication that 
the photoionization cross section of an endohedral atom 
may be rather insensitive to the degree of diffuseness of 
the confining potential boarders, in reasonably large lim- 
its. As a note on an independent issue, one can see that 
the photoionization cross section of free hydrogen dif- 
fers from that of HQCgo by the presence of a very well 
developed maximum at about 25 eV and another less de- 
veloped maximum at about 43 eV in the latter. These 
are called confinement resonances, i.e., resonances that 
occur due to a constructive interference of the outgoing 
photoelectron wave and photoelectron waves scattered off 
the Ceo cage [l|, H, Q (and references therein) . Confine- 
ment resonances have only recently been experimentally 
proven as reality ^3|]. 

In the next, we investigate whether the above findings 
are valid for photoionization of inner-shells and valence- 
shells of endohedral multielectron atoms A@Gqo as well. 
To date, only photoionization of endohedral Xe@C^Q 
has been reliably measured experimentally Q. There- 
fore, we choose photoionization of XeOCeo as the case 
study to learn how discontinuity of J/swpC?") and diffuse- 
ness of L''dp('') niay affect photoionization cross sections 
(T„£(a;) of, and photoelectron angular-asymmetry param- 
eters l3ni{io) from, the XeQCgo inner 4d^° and valence 
5p^ subshells. Note, Pne{uj) depends on phase shifts Si±i 
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FIG. 3: (Color online) RPAE calculated data for 0-4^(0;) and 
/?4(i(a;) of XeQCeo obtained with the use of the square- well- 
potential ?7swp('") (dashed lines) and diffuse potential Unp{r) 
with 77 = 0.1 (solid lines), as well as corresponding IfF cal- 
culated data obtained with the use of UDp{r) with ?7 = 0.1 
(dotted lines), as marked. Also depicted are corresponding 
RPAE calculated data for free Xe (dashed-dotted lines) and 
experimental data for a4d{uj) of Xe@C^o 01 ■ Experimental 
data were shifted by 4.2 eV towards higher photon energies 
as well as multiplied by 10 to compare with theory (see dis- 
cussion in Appendix). 
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FIG. 4: (Color online) RPAE calculated data for XeOCeo 5p 
photoionization cross section, (T5p{uj), and dipole photoelec- 
tron angular-asymmetry parameter, I3^p{u)), calculated with 
the use of the square- well-potential !7swp('') (dashed lines) 
and diffuse potential U-Dp{r) with 77 = 0.1 (solid lines). Also 
depicted are RPAE calculated data for a^piuj) and I3^p{u}) for 
free Xe (dashed-dotted lines). 



III. CONCLUSION 



of photoionization matrix elements. Hence, correspond- 
ing calculated data will provide, although implicitly, the 
information on sensitivity of phase shifts 5i±i to difFuse- 
ness of a confining potential as well. RPAE calculated 
data for (T4d(u;) and /34d(a;), as well as for azp{Lu) and 
P^p{oj), are depicted in Fig. [3] and Fig. SI respectively. 

One can see from Figs. [3]and|4] (compare dashed curves 
marked as "/7swp" with solid curves marked as "J7dp, 
r] = 0.1" ) that, similar to the HQCeo case, there are little- 
to-no differences both in calculated CT„f (a;)'s and /3„f (lli)'s 
obtained with accounting for a square- well-potential con- 
finement versus those obtained with a diffuse confine- 
ment potential even with as large degree of diffuseness as 
r] — 0.1. Thus, photoionization spectra both of a single- 
electron and multielectron endohedral atom AQCgo are 
equally insensitive to discontinuity in J/svvp or diffuse- 
ness of C/dp; the implication is that the findings carry 
the general meaning. 



To summarize, it has been proven in this paper that 
discontinuity of a square-well-potential does not lead to 
any artifacts, whether quantitative or qualitative, in pho- 
toionization spectra of endohedral atoms. The square- 
well-potential model of Cgo is as good as its modeling 
by a diffuse potential. Both of these potentials lead to 
practically identical calculated data for photoionization 
spectra. The latter appear to be largely insensitive to the 
degree of diffuseness. Hence, either of the two potentials 
is equally suitable for mimicking the Cgo cage. In par- 
ticular, switching from a confining square-well-potential 
to a more realistic diffuse potential, if wanted, will not 
require a re-study of a rich variety of important results 
on A@Gqo ionization phenomena obtained to date on the 
basis of a confining square- well-potential concept. 

Finally (see Appendix), since the values of square- 
well-potential parameters i?o — 6.01, A — 1.25, and 
Uq — 0.422 a.u., empirically deduced in Ref. [2^ and 
exploited in the present work, lead to a reasonable agree- 
ment with experiment (see Fig. ^ , we suggest the users 
of the square-well-potential modeling of A@Cqq to uti- 
lize these updated parameters in their future work on 
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the subject. 
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Appendix 

To finahze the paper, we now briefly comment on par- 
ticular features and structures in the calculated spectra 
of XeCSCgo depicted in Figs. |3] and IH since they are im- 
portant stand-alone qualities. 

1. The 4d photoionization of Xe@C60 

To start with, the three strong oscillations in a4d{^) 
of Xe@C6o about fT4d(w) of free Xe are confinement reso- 
nances (see the definition of confinement resonances given 
above in this paper) . They were previously theoretically 
studied and detailed in Refs. [1, H, El [ll ill and ex- 
perimentally in Ref. [3] . Results of the measurement are 
by about the factor of 10 smaller than the cited theoret- 
ical data. This is because the measured channel of the 
Xe@Cgo id photoionization might [1] account for only 
about 10% of the Xe id oscillator. For this reason the 
experimental data depicted in Fig. |3] of the present pa- 
per are multiplied by 10, to ease the comparison with 
theory. One can see that calculated results for a4d{uj) of 
XeQCeo obtained in the present work with the square- 
well-parameters Rq = 6.01, A = 1.25, and Uq = 0.422 
a.u. are in a reasonable agreement with main structures 
of the experimental spectrum, to a good approximation. 
Interesting, HF calculated data for (T4d{uj) of Xe@C6o 
seem to be even in a better, than RPAE data, agreement 
with experiment beyond the photon energy of about 100 
eV. Note, the comparison of our calculated data for neu- 
tral XeHCeo with experimental data for charged Xe@C(iQ 
is appropriate because P, [13] charging the Ceo shell 
positively does nothing to the photoionization cross sec- 
tion (as a function of photon energy) except to increase 
the threshold energy. 

Results for depicted in Fig. [3]/34(i(w) of XeQCgo illus- 
trate the impacts of the confining square- well and diffuse 
potentials (when compared with free Xe data) as well as 
electron correlation (when RPAE is compared with HF 
data) on the P4d{uj) spectrum. They are self-explanatory. 
We omit commenting on them. 

To conclude the discussion on the XeQCgo id pho- 
toionization we note that very recently [Tgf new theo- 
retical data on the XedCgo id photoionization were re- 
ported. Although being obtained on the basis of exactly 



the same approximations as in the works of three other 
independent theoretical grou£S [H [H HI] [i.e., RPAE 
(or relativistic RPAE in Ref. [l^) & square- well-potential 
with identical parameters) they differ strongly from re- 
sults of all the three cited theoretical groups (the latter 
are in the mutual agreement with each other). We con- 
clude that something was done wrong in work 19], most 
likely because of a very peculiar way those authors chose 
for solving HF equations in the presence of a square- well 
potential confinement. 



2. The 5p photoionization of XeQCeo 

Commenting on cr5p(w) of XeQCeo, Fig. I^lja), we note 
that a very strong maximum at the 5p threshold as well 
as a much less developed oscillation at about 30 eV carry 
the confinement resonance origin. Interesting, they are 
absent in /35p(a;). 

As for a5p{uj) beyond of a 100 eV energy, it domi- 
nates by the three strong resonances. Their positions 
approximately match the positions of confinement res- 
onances in a4d(uj), Fig. [31 Exactly the same situation 
has recently [29| been predicted and explained for the 
XeQCgo 5p generalized oscillator strength in a fast elec- 
tron impact ionization of XeQCgo. As in Ref. [1^, the 
present study reveals that the resonances in a5p{oj) be- 
yond 100 eV are induced by the confinement resonances 
in the id photoionization channels, via interchannel cou- 
pling. Thus, the resonances in a^plcu) beyond 100 eV are 
of the correlation confinement resonance nature. Corre- 
lation confinement resonances in AQCgo photoionization 
[E [3 [IE HqI are resonances that emerge in the pho- 
toionization of an outer subshell due to interference of 
transitions from the outer subshell {5p in our case) with 
confinement resonances emerging in inner shell transi- 
tions (id transitions in our case), via interchannel cou- 
pling. These correlation confinement resonances in the 
XedCeo 5p spectra are strong because of, as in free Xe, 
a strong intershell correlation coupling between id on 5p 
transitions. 

Furthermore, one can see from Fig. IDJb) that corre- 
lation confinement resonances show up prominently in 
/35p(a;) of XeQCgo as well [resonances at about 85 and 
100 eV in /35p(a;)]. This is another interesting prediction 
of the present study. 

Finally, we note, that we omitted the Cgo dynamical 
polarization impact (lol . [3l| (and references therein) on 
the XeQCgo 5p ionization which is known to be signifi- 
cant up to about 70 eV of photon energy. The omission 
is justified because the primarily aim of the present pa- 
per is to study differences between effects of the confin- 
ing square- well-potential and diffuse potential on A@Cgq 
photoionization, and the Cgo dynamical polarization can- 
not alter said differences. 
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